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Abstract 

In this paper we give all the details of the calculation that we presented 
in our previous paper [arXiv:0712l 3522, concerning the four-loop anomalous 
dimension of the Konishi descendant tr(0Z0Z — (jxpZZ) in the SU (2) sector 
of the Af = 4 planar SYM theory. We explicitly consider all the wrap- 
ping diagrams that we compute using an A/" = 1 superspace approach and 
Gegenbauer polynomial x-space techniques. 
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1 Introduction 



Recently we have presented [1] the calculation of the anomalous dimension of the com- 
posite operator ii{(j)Z(j)Z — (pcpZZ) at four loops in the Af — 4 planar supersymmetric 
Yang- Mills theory, where and Z are two complex scalar fields of the theory. 

The importance of this calculation resides in the fact that this is the simplest case 
in which the so-called wrapping effects are present: at a given perturbative order K in 



where g — gyu is the Yang-Mills coupling constant, the range of the interactions between 
adjacent fields grows with the perturbative order as i^-l- 1. For an operator with L fields 
we should expect new effects when the range exceeds L, that is at order g"^^. It turns 
out that L = 4 is the simplest case in which these effects can be seen (operators with 
L — 2 and L — 3 are automatically BPS). 

Understanding these wrapping effects is a crucial step towards a comparison between 
the spectrum of the anomalous dimensions of gauge invariant operators of planar A/" = 4 
SYM and the spectrum of strings on AdS5 x S^, as predicted by the AdS/CFT conjecture 



Indeed, in the case of long operators, the comparison between the spectra of the two 
theories is now at hand, thanks to the recent progress in understanding their integrability 
properties. On the gauge theory side it was crucial the realization [3,4] that the planar 
one-loop dilatation operator of A/" = 4 SYM maps into the Hamiltonian of an integrable 
spin chain. This spin chain picture was extended to higher orders in perturbation 
theory [5, 6] and thanks to the integrability properties of the system the computation 
of the anomalous dimensions was then reduced to finding solutions of the associated 
Bethe equations [7] . This program cannot be applied to short operators since the Bethe 
ansatz is only asymptotic and wrapping effects appear at order g"^^: the interaction is 
no longer localized in some limited region along the state and asymptotic states cannot 
be defined. 

Clearly from the perspective of integrable models including wrapping effects is very 
important and several papers have addressed this issue. In [8] the properties of wrap- 
ping interactions have been analyzed in terms of Feynman diagrams. In [9] wrapping 
interactions in the BMN matrix model have been discussed. In [10, 11] it was proposed 
to study wrapping interactions via the thermodynamic Bethe ansatz, while in [12] it was 
assumed that wrapping effects might be described by the Hubbard model. In [13] the 
form of wrapping contributions at four loops was conjectured by using restrictions from 
the BFKL equation [14-16]. In [17] wrapping effects in some simple spin chain models 
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were studied. On the string theory side, finite size contributions have been analyzed in 
several recent papers [18-25]. 

In order to test these proposals it is crucial to perform an explicit field theoretical 
computation of an anomalous dimension at an order in which wrapping contributions 
become relevant. As mentioned above the simplest such a case is the four- loop anomalous 
dimension of the Konishi descendant ti {(f)Z(f)Z—(j)(f)ZZ) in the SU{2) sector of the theory. 

In this paper we describe in detail the computation we have presented in [1], with 
the aim of supplying the interested reader with all the details needed to reproduce the 
calculation. 

The complete four-loop Feynman graph calculation is very complicated. However one 
can take advantage of the approach based on the asymptotic dilatation operator and 
avoid the explicit computation of some of the most difficult classes of graphs. Indeed 
the contribution from all non-wrapping graphs, corresponding to interactions of range 
from one to four, is included in the four-loop asymptotic dilatation operator together 
with the interactions of range five. 

Then one is left with two operations to be done: the subtraction of the range- 
five contributions from the dilatation operator, and the computation of the wrapping 
graphs. These two operations are drastically simplified by means of A/" = 1 superspace 
techniques. The integrals that we produce after D-algebra manipulations on the su- 
pergraphs are computed by the Gegenbauer polynomial x-space technique (GPXT) [26] 
and by the method of uniqueness [27]. 

The anomalous dimension of a composite operator O is extracted from the - pole of 
the graphs contributing to its renormalization in a standard way. In the special case in 
which the operator is multiplicatively renormalized the anomalous dimension is given 
by 

7(0) = lim Eg-^\ogZai9,e) , (1.2) 
e-*o I dg J 

where 

Cren = C'bare • (1-3) 

In the general case of operator mixing, the last equation should be rewritten in matrix 
form and the anomalous dimensions are given by the eigenvalues of the matrix appearing 
on the right hand side of (11.21) . In our case we have two mixing operators but it turns out 
that the combination tr(0Z0Z — (fxpZZ) is an exact eigenvector of the renormalization 
matrix to all orders, so it is multiplicatively renormalized. 

The most remarkable outcome of our calculation is that at four-loop order wrapping 
interactions give rise to contributions proportional to ({5) increasing the level of tran- 
scendentality of the anomalous dimensions. Our result rules out the previous conjectures 
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on how to treat wrapping contributions proposed in [12,13,28]. 

Another computation of the same anomalous dimension appeared later on in [29]. 
The results presented there do not match our findings, but the authors did not attempt 
a comparison with the computations we had reported. The result of [29] contains a 
trivial error in the final step, however even after correcting it, the result is different from 
ours. Finding the origin of the discrepancy is not an easy task. The authors of [29] 
used exactly the same strategy as ours but their calculation of the wrapping diagrams is 
performed using a component field approach. Thus the number of diagrams they have to 
compute is huge, supersymmetry is not maintained explicitly, renormalization is subtle 
and partial comparison with our classes of diagrams is almost impossible. Needless to 
say that keeping under control such a large number of contributions is very difficult. On 
the contrary using superfields the calculation is much simpler since several cancellations 
are automatic, and the number of relevant diagrams is significantly reduced. In any 
event we hope that all the detailed explanations we are going to present in this paper 
will help to clarify the situation. 

Our paper is organized as follows. In Section [2] we recall the basic ingredients of the 
Af = 1 superspace approach, and we discuss an important result on the displacement of 
spinorial derivatives on a supergraph that allows us to drastically reduce the number of 
potential contributions. In Section [3] we describe how to subtract the range-five terms 
from the asymptotic Hamiltonian and we build the effective Hamiltonian containing all 
the contributions from range one up to range four. Thanks to a cancellation described in 
Appendix \^ the range- five subtraction procedure can be performed without computing 
any Feynman diagram. However in order to provide a further check of our procedure 
we have computed all range-five graphs explicitly. This is reported in Appendix [Bl 
Subsection 14.11 contains the calculation of the wrapping supergraphs. The full list of 
them is reported in Appendix [Cl with a description of all the cancellations occurring 
among them and with all the integrals produced after completion of the D-algebra. The 
GPXT method used to compute the momentum integrals is summarized in Subsection 
14.21 and Appendix [Dl Using the results obtained so far in Subsection 14.31 we are able 
to give an expression for the four-loop wrapping dilatation operator. In Section [5] we 
collect our results and we give the final answer for the Konishi anomalous dimension at 
four loops. 

2 Af = I superspace approach 

We are going to describe the calculation of the anomalous dimension of the composite 
operator tr(0Z0Z — (jxpZZ) at four loops in the planar limit. In order to streamline the 
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calculation we decided to work using a A/" = 1 superspace formalism. In this section we 
review the most relevant features of this approach needed to perform our calculation. 

The action of A/" = 4 SYM can be written using Af = 1 superspace in terms of one 
real vector superfield V and three chiral superfields 0* (we follow the notations and 
conventions of [30]) 



S = j d^x d^5 tr (e-^'^ e^'^ '^O + ^ J ^^^^ tr {WW^ 
+ 1^1 d^xd^ee'^' tr(0, [0,,0fc])+h.c. , 



(2.1) 



where = iD^ (e"^^ e»^) , and 1/ = VT'', 0* = 0^T'^, i = 1, 2, 3, T" being matrices 
satisfying the SU{N) algebra 

[Ta,n]=tUcTc (2.2) 

and being normalized as 

tr(r,r,) = 6at . (2.3) 

We will usually denote the three chiral superfields 0* as {4>,ip, Z). 

We want to study the renormalization of the composite chiral operators 

Ci = tT{(j)Z(f)Z) , O2 = tT{(f)(j)ZZ) , (2.4) 

which mix under renormalization. To this end, the first step is the construction of all 
the supergraphs containing one operator insertion. The Feynman rules can be derived 
from the action fl2.ip . In momentum space we have the superfield propagators 

while the vertices are read-off directly from the interaction terms in (12. ip . with additional 
l)^, D"^ factors for chiral, antichiral lines respectively. The ones that we need are the 
following 

a \ (2.6) 

Vertices containing three or more vector superfields V do not enter our calculation. 

The general idea is to consider a given supergraph and to perform the Z)-algebra in 
order to reduce it to a standard graph. In doing so one usually produces a large number 
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of contributions, but most of them are not relevant for the computation of the anomalous 
dimension of the composite operator since they give rise to finite integrals. Thus it is 
very useful to know and to exploit D-algebra manipulations that can be made on the 
supergraph and that are legal only up to finite terms. In the next subsection we are 
going to describe a property of the supergraphs that was very important in computing 
our diagrams. 

2.1 A useful proof 

We consider diagrams with one operator insertion, built using only chiral, antichiral, 
single-vector and double-vector vertices. We arc interested in isolating those contribu- 
tions which give rise to divergent integrals, i.e. the ones which will contribute to the 
anomalous dimension. Here we want to prove that in order to produce divergent con- 
tributions we have to perform and complete the D-algebra in such a way that no spinor 
derivative is moved out of the diagram onto the external lines, except for derivatives 
on scalar propagators that do not belong to any loop from the start. Let us define the 
following quantities: 

• Vc — number of chiral vertices 

• Va — number of antichiral vertices 

• Vy^ — number of single- vector vertices 

• Vy = number of double-vector vertices 

• Pc = number of scalar propagators belonging to at least one loop 

• Py = number of vector propagators 

• p = total number of propagators belonging to at least one loop 

• E = number of scalar propagators not belonging to any loop. 

• Nj( — number of loops 

• Nr,, Nf)-. number of spinor derivatives 

In all our diagrams we have Va — Vq and Njj — N^,, and E is equal to the number of 
chiral vertices not belonging to any loop. 

Moreover the number of outgoing external lines equals the number of fields in the com- 
posite operator. From each chiral or antichiral vertex three scalar lines start. From each 
single- vector vertex, two scalar and one vector propagators start. Prom each double- 
vector vertex, two scalar and two vector lines start. 
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Given these relations, we can write 



3{Vc + Va)+2{V^'^ + V;f>) 



(2.7) 



At each chiral vertex, two out of the three scalar lines have a D^. 

At each antichiral vertex, two out of the three scalar lines have a D"^. 

At each single- vector or double- vector vertex, one of the scalar lines has a D'^, while the 

other one has a, D"^. So we have 



Then we can write 

p =Pc +Pv 



a) 



(2)^ 



No = AVa + 2{V^^^ + V, 



K2)^ 



Vc + Va + Nd + v. 



(1) 



V 



2K 



(2) 



V 



-E 



(2. 



(2.9) 



In order to obtain a logarithmic divergent diagram, the £)-algebra must produce a 
number of standard derivatives equal to 



2p - 4iV£ 



(2.10) 



We need a D (and a D) to create each standard derivative. Furthermore, we need 2A^^ 
extra D's (and D's) to complete the superspace d^^ integration for the loops. So 
after performing the Z^-algebra, a momentum integral with a surface divergence can be 
found only if at least 

2p-2Ne (2.11) 

out of all the £)'s are not brought outside the diagram on the external lines. 

A can always be moved onto each scalar propagator which is not part of a loop: 
either it is already there at the beginning of the D-algebra, or it can be moved there 
through an integration by parts at the internal vertex at which the external propagator 
is linked. So the effective number of L>'s available in the loops is 



Nd-2E 



(2.12) 



During the D-algebra manipulations we will be allowed to move some D^s outside the 
diagram only if their number at the start exceeds the number we need in order to produce 
a divergent integral. So we should have 



Nd-2E >2p- 2Ne 



Vc + Va + Nd + Vj'^ + 2V^') 



2E - 2N, 



(2.13) 
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that is 



Np> - 



(2.14) 



This condition is never satisfied by the diagrams in the class we are studying, as can be 
seen using Euler's formula for connected graphs 



V -P + F 



(2.15) 



where V denotes the number of vertices, P is the number of propagators and F is the 

number of faces. 

In our case (remember that the composite operator behaves as an additional vertex) we 
have 



which leads to 



No 



1 

2 L 



Vc + Va + + 2\4') 



F = + 1 , (2.16) 



(2.17) 



We conclude that in order to produce a divergent contribution no spinor derivatives can 
be brought outside the diagram. 

All our D-algebra manipulations were performed making use of this property: only 
the contributions in which all spinor derivatives stay in the loops were kept. Clearly 
this rule allowed to discard a huge number of irrelevant contributions. 



3 Subtraction of range-five diagrams from the asymp- 
totic Hamiltonian 

In this section we are going to describe how to compute the contributions from the 
diagrams with range from one to four with no wrapping interactions. In order to avoid 
their explicit, very complicated Feynman graph computation we take advantage of the 
asymptotic four-loop dilatation operator D4, whose expression was given in the SU{2) 
sector by [31]: 
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(3.1) 



Di= - (560 + 4/?){} 

+ (1072 + 12/5 + 8e3a){l} 

- (84 + 6/? + 4e3a){l, 3} - 4{1, 4} - (302 + 4/3 + 8e3a)({l, 2} + {2, 1}) 
+ (4/? + 4e3a + 2ze3e - 42e3d){l, 3, 2} + {A/3 + 4e3a - 2tesc + 4ze3d){2, 1, 3} 
+ (4 - 2ze3c)({l, 2, 4} + {1, 4, 3}) + (4 + 2^e3c)({l, 3, 4} + {2, 1, 4}) 
+ (96 + 4e3a)({l,2,3} + {3,2,l}) 

- (12 + 2/? + 4e3a){2, 1, 3, 2} + (18 + 4e3a)({l, 3, 2, 4} + {2, 1, 4, 3}) 

- (8 + 2e3a + 2ze3b)({l, 2, 4, 3} + {1, 4, 3, 2}) 

- (8 + 2e3a - 2ze3b) ({2, 1, 3, 4} + {3, 2, 1, 4}) 
-10({l,2,3,4} + {4, 3,2,1}) , 

where e^a, ^3b, ^sc and 63^ parameterize the invariance of the Hamiltonian under similarity 
transformations. In a perturbative approach their values depend on the choice of the 
renormalization scheme. The parameter (3 = 4(^(3) comes from the dressing phase 
[28,31-34]. 

The permutation structures appearing in fl3.ll) are defined as 

L-l 

{^1, . . . , a„} = Pai+r ai+r+l ■ ■ ■ Pa„+r a„+r+l ; (3-2) 

r=0 

where Paa+i, with the cychc identification Paa+i — Pa+L a+L+i, permutes the fiavours 
of the ath and (a + l)th site when acting on a cyclic state of length L. The number of 
nearest neighbours, interacting in the fiavour space, is extracted from the list of integers 
tti, . . . ,an in the permutation structure as 

K = 2 + max — min . (3.3) 

ai...an ai...an 

It is limited from above by the order K in the perturbation expansion in powers A^^ as 
K, < K + 1. The asymptotic regime in which D4 can be trustfully applied is determined 
by the relation L > K. Some rules valid in the asymptotic case for the manipulation of 
these structures can be found in [35]. 

The expression f l3.ip contains contributions which describe the permutations among 
K = 5 neighbouring legs. Therefore it is correct only if applied to a state in the asymp- 
totic sense, i.e. the number of sites in the state has to be five or more. We can correct 



8 



it for the application on a length four state: the contributions from all Feynman dia- 
grams which describe the interactions of five neighbouring legs have to be replaced by 
the contributions from all four-loop wrapping interactions. 

Before doing this, it is advantageous to change the basis from the permutation struc- 
tures (I3.2p to a set of functions which are directly related to the chiral structures of the 
Feynman supergraphs. At four loops we are lead to define the following chiral functions 



(3.4) 



X(a, b, c, d) = {} - 4{1} + {a, b} + {a, c} + {a, d} + {b, c} + {b, d} + {c, d} 
— {a, b, c} — {a, b, d} — {a, c, d} — {b, c, d} + {a, b, c, d} , 
X(a, b, c) = -{} + 3{1} - {a, 6} - {a, c} - {6, c} + {a, 6, c} , 
Xia,b) = {}-2{l} + {a,b} , 
= -{} + {!}, 

x() = {}. 

The relevance of these functions stems from the fact that any Feynman supergraph 
produces a permutation structure given by one of them. Moreover we call them chiral 
because their form depends only on the structure of the chiral interactions of the dia- 
gram, while they are completely insensitive to the presence of vector lines. 
The form of these functions can be easily derived: xO is just the identity, corresponding 
to a Feynman supergraph without chiral vertices. The first non-trivial function, x(l) 
corresponds to the diagram given in Figure [H 



+ 



z 



V 

i 

A 



Figure 1: Building block 



This diagram has two vertices V3 and V3 (see (12.61) ) connected by a {ipTp) superfield 
propagator, which is killed by the D-algebra. It is a sort of effective four-vertex, very 
similar to the one describing the scalar interaction in a component approach. This dia- 
gram is a "building block" for the construction of more complicated ones. 
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Indeed it is easy to see that every four-loop supergraph contributing to the one-point 
function of a chiral operator in the SU{2) sector, hke our operators fl2.4p . are built by 
gluing together up to four of the building blocks given in Figure [TJ The form of the 
chiral functions in (13 ■4p then follows consequently. 

Let us notice that the number of arguments n, with n = 1, ... ,4 at four loops, corre- 
sponds to the number of building blocks in our underlying Feynman graphs. The fact 
that the coefficients in each chiral structure sum up to zero guarantees that the two 
length L BPS operators ti^Z^) and tr(0Z'^~^) are protected. 

In order to rewrite in (13. ip in terms of the chiral functions, we have to express 
the permutation structures (13. 2p in terms of the chiral functions (13. 4p . We have 



Making use of these relations we reexpress the four-loop dilatation operator (13. ip in 
terms of the chiral functions and obtain 



D,= + 200x(l) - 150[x(l, 2) + x(2, 1)] + 8(10 + e3a)x(l, 3) - 4^(1, 4) 
+ 60[x(l,2,3)+x(3,2,l)] 

+ (8 + 2/? + 4e3a - 4ze36 + 2^63^ - i^e3d)x{l, 3, 2) 
+ {8 + 2(3 + 4e3a + 42636 - 22e3c + 42e3d)x(2, 1, 3) 

- (4 + Ate,, + 2ze3e) 2, 4) + x(l, 4, 3)] 

- (4 - 4ze3fe - 2tesc) 3, 4) + x(2, 1, 4)] (3.6) 
-(12 + 2/? + 4e3a)x(2,l,3,2) 

+ (18 + 4e3.)[x(l, 3,2, 4) + x(2, 1,4,3)] 

- (8 + 2e3a + 2test) 2, 4, 3) + 4, 3, 2)] 

- (8 + 2e3a - 2^636) [X(2, 1,3,4) + x(3, 2, 1, 4)] 
-10[x(l,2,3,4) + x(4,3,2,l)] . 



This representation is particularly useful since the chiral functions x and not the permu- 
tation structures are directly related to the underlying Feynman supergraph. However 



{a,b,c, 



d} 



X{a, b, c, d) + x{a, b, c) + x{a, b, d) + x{a, c, d) + x{b, c, d) 
+ x(a, b) + x(a, c) + d) + xip, c) + xip, d) + x(c, d) 
+ 4x(l) + x() , 

X(a, b, c) + x(a, b) + x(a, c) + xip, c) + 3^(1) + xO , 
X(a,6)+2x(l) + x() , 

x(i) + x() , 



{a, 6, c} 
{a, 6} 
{1} 
{} 



(3.5) 
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there is still an ambiguity: two supergraphs which differ only by the arrangement of 
the flavour-neutral vector superfield lines contribute to the coefficient in front of the 
same chiral function. Therefore the chiral functions do not allow in general to extract 
information about the range of the interaction of the underlying Feynman supergraph. 

In particular range-five diagrams in which the first or fifth line is connected to the 
rest of the diagram only via vector line(s) appear with the same chiral function as the 
corresponding range-four diagram in which the corresponding first or fifth line is a non- 
interacting spectator. We call those range-five diagrams which have a chiral function 
which does not capture the range of their interaction non-maximal range-five diagrams. 
They also have to be subtracted from D4 and in general one should compute explicitly all 
these Feynman graphs. Luckily it turns out that there is no need to do that. In Appendix 
|X]we show that the non-maximal lange-fLve interactions do not contribute to -D4. This 
means that we can determine the subtraction directly from D4. This operation becomes 
particularly simple using the representation of in terms of the chiral functions x- 
We only have to neglect all the terms with chiral functions x{C'iy • • • ; (^n) which describe 
a range-five interaction, i.e. which fulfill max — min = 3. Hence we have to subtract 

ai...a„ ai...a„ 

from (13.61) 

6D, = - 10[x(l, 2, 3, 4) + x(4, 3, 2, 1)] + (18 + 4e,a) 3, 2, 4) + x(2, 1, 4, 3)] 

- (8 + 2e3„ + 2^635) 2, 4, 3) + x(l, 4, 3, 2)] 

- (8 + 2e3„ - 2tesh) [x(2, 1, 3, 4) + x(3, 2, 1, 4)] 

- (4 + 4ze3b + 2ze3e) 2, 4) + x(l, 4, 3)] 

- (4 - Mesb - 2zesc) 3, 4) + x(2, 1, 4)] 
-4x(l,4) . 

We should stress that due to the decomposition in fl3.4p . the subtraction is not 
only a simple neglection of all range-five permutation structures as attempted in [9] for 
the BMN matrix model. One has to modify also the coefficients of the permutation 
structures of lower range in (13.11) . 

In order to have an independent check of the whole procedure, we have found the 
above constructed range-five contributions (13.71) also through an explicit Feynman graph 
computation. We give in Appendix [B] all the details of the calculation. As a byproduct 
we have found in this way explicit values for the unphysical coefficients e^a, est, esc- In 
the MS-scheme we have used, they are given by 

4 4 

eSa = -4 , 636 = , e3c = I- . (3.8) 
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Thanks to the absence of non-maximal range-five interactions the explicit Feynman 
graph computation was not necessary. It was performed just as a further check on the 
range- five part of D4. 

We emphasize that a straightforward modification of the proof in Appendix |X] leads 
to the cancellation of this type of contributions also at higher loop orders. This feature 
is quite interesting since it allows to determine the necessary subtractions in the general 
case of Dk applied to a length L = K operator, with no need of Feynman diagram 
computations. 

Finally after subtraction of the range-five contributions the four-loop dilatation op- 
erator becomes 



We can now apply the subtracted dilatation operator to the states with L = 4 sites. In 
the SU (2) subsector there exist two composite operators of length L = 4 and with two 
'impurities' which mix under renormalization. The corresponding states are given by 



X(1)^M, ;^(l,2)^-M, x(l,3)^-2M, ;^(2, 1) -M , 
X(1,2,3)^M, x(3,2,l)->M, x(2,l,3)-2M, 3, 2) 2M , 



/^sub ^ _ ^ + 200x(l) - 150[x(l, 2) + x(2, 1)] + 8(10 + e3.)x(l, 3) 



+ 60[x(l,2,3) + x(3,2,l)] 

+ {8 + 2(3 + 4e3a - 4^636 + 2ze3c - 4^e3d)x(l, 3, 2) 
+ (8 + 2/5 + 4e3a + 4ze3b - 2ie-ic + 4^e3,)x(2, 1, 3) 
-(12 + 2/? + 4e3a)x(2,l,3,2) . 



(3.9) 



Oi = tr(0Z0Z) , O2 = tT{^4>ZZ) . (3.10) 

Defining a two-dimensional vector O = {Oi, 02)^, the application of the chiral structures 
to O is captured by the replacements 



X(2,l,3 



2) ^ -2M 



(3.11) 



where the common mixing matrix is given by 




(3.12) 



The eigenvectors and eigenvalues of M are given by 




(3.13) 
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These are directly the eigenvectors of in fl3.9l) . since in the L = 4 operator basis 
(I3.10p it is proportional to M. With /? = 4^(3) and e^a from fl3.8l) it assumes the form 

Df^ ^ 4(121 + 12C(3))M . (3.14) 

In the next section we move our attention to the calculation of wrapping diagrams. 

4 Contributions from wrapping diagrams 

In this section we are going to face the computation of the wrapping contributions. 
This is the part of the analysis that necessarily requires the explicit computation of 
Feynman supergraphs. The first step will be to build all wrapping supergraphs, and for 
this we rely on the results of [8] , where a systematic Feynman-diagrammatic analysis of 
wrapping interactions was performed. 

In the present paper we are using Af = 1 superspace techniques, so we need to adapt 
that approach for the construction of all the wrapping supergraphs contributing to the 
renormalization of the chiral operators fl3.10p . 

The second step will be to compute the four-loop momentum integrals using the 
GPXT. 

Finally we collect these results in a wrapping four-loop dilatation operator. 
4.1 Wrapping supergraphs 

First of all we have a set of wrapping graphs with only chiral interactions, which are 
shown in Figure [21 




Figure 2: Wrapping diagrams with only chiral interactions 



With the identification of the first and the fifth lines, we can still use the defini- 

13 



tion fl3.2p for the permutation operators. Then the chiral structures will be 



Wax ^ X(2,4,l,3) ^ -4M 

Wa2 ^ X(4,l,2,3) ^ -M (4.1) 
Wa3 ^ X(4,3,l,2) ^ -2M 

The reflection of Wa2-, which is not symmetric, will have the structure x(l5 4,3,2) — > 
-M. 

Beware that, unlike in the asymptotic case, here the exchange 1 4 is not allowed, 
since lines 1 and 4 are now neighbours. Furthermore, there is an ambiguity in the choice 
of the basis. Shifts of all arguments by an integer generate four different chiral functions 
for Wai and Waz and two each for Wa2 and for its reflection. A cyclic rotation of the 
external legs of Wa^^ generates an additional chiral function. The elements in each of 
the four classes are equivalent when applied to a length-four state. From each class we 
have chosen a single representative. 

The results after completion of the D-algebra for each completely chiral supergraph 
are given by 

Wai ^ {g^N)\h/2) x(2, 4, 1, 3) ^ -2{g''NfhM , (4.2) 
Wa2 ^ (/iV)4/2[x(l,4,3,2) + x(4,l,2,3)]^-2(/iV)4j2M, (4.3) 
Waz ^ (/iV)%x(4,3,l,2)^-2((7'iV)%M, (4.4) 

where /j denote the momentum integrals listed in Table IC.8[ 
Then the total contribution from this class reads 

Wa.. ^ -2{g^N)\h + I2 + h)M . (4.5) 

Graphs with chiral structure x(l)2,3) are shown in Figure ICTTI while the results 
of the D-algebra for each diagram are summarized schematically in Table IC.ll here 
cancellations between pairs of diagrams, as well as graphs which do not produce any 
divergence are explicitly indicated. Diagrams which do not cancel and therefore must 
be fully computed are marked with a *. For these relevant graphs we also give explicitly 
the value of the relative symmetry factor. Graphs with chiral structure x(3, 2, 1), which 
is the reflection of x(l, 2, 3), contribute exactly in the same manner. We give a summary 
of the full results for all the structures in Table [Q7l where the explicit expressions (13.111) 
and fl4.ll) for the chiral structures in the length-four subsector have been used. 
So the total contribution from these chiral structures is given by 

Wb*. ^ -A{g^N)\h - h)M . (4.6) 
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We proceed in the same way for the other structures. Diagrams associated to 
x(l,3,2) are hsted in Figure [ni2l while Table \C2\ contains the results from D-algebra. 

J2Wc**^-^9'NYI,M . (4.7) 

There is one more structure which needs only a vector to be completed to four loops, 
x(2, 1,3). The corresponding diagrams and the D-algebra results can be read from 
Figure [U73] and Table [Cl3l The divergent contribution i^ 

J2Wd»^~^9'NYIsM . (4.8) 

Next we turn to chiral structures of shorter range. Two of them require two vectors 
to be completed to four loops: the first one is x(2, 1), for which diagrams and Z)-algebra 
results are given in Figure IC.4I and Table IC.4[ Adding the analogous result from graphs 
with structure x(l,2) we obtain 

J2 We*. ^ -A{g^N)\lr - h)M . (4.9) 

The second structure to be completed by two vectors is x(l)3) it is described in 
Figure IC.5I and Table IC.51 from which we get 

Wf** ^ -4{g^N)\ho + In - 2/i2)M . (4.10) 

The last chiral structure to be considered is x(l), which is the only one requiring three 
vectors. The diagrams for this structure are presented in Figure IC.6I Using the results 
summarized in Table IC.6I we obtain 

^ Wg.. 2{g^N)\h - h)M . (4.11) 

So up to this point we have obtained the total contribution from each relevant chiral 
structure. Now we have to substitute the values of the momentum integrals In order 
to obtain the full wrapping contribution to the dilatation operator on the length-four 
subsector, we have to take the coefficient of the ^ pole of these expressions and multiply 
them by —8 according to eq. fll.2p . In the next subsection we describe the computation 
of these momentum integrals by means of the Gegenbauer polynomial x-space technique. 

^In the first version of tlic paper tliere is a symmetry factor of 2 missing in this expression. We have 
found this factor after a revision of our calculation which was suggested by a mismatch between the 
rational part of our previously found Konishi anomalous dimension and the one computed in [36] . 
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4.2 Gegenbauer polynomial x-space technique 

The integrals with non-trivial numerators which cannot be rephrased in terms of stan- 
dard scalar integrals, are calculated with the Gegenbauer polynomial x-space technique 
(GPXT) [26], and are independently checked with the method of uniqueness [27] and 
with the help of MINCER [37], a computer program to compute 3-loop integrals. We have 
furthermore used the GPXT and MINCER to reproduce the known results for the scalar 
integrals. 

The GPXT is based on the fact that in x-space integrals all propagators depend on 
the difference of only two coordinates and hence can be expanded in a series involving 
the Gegenbauer polynomials. The non-trivial numerators given by scalar products of 
several momenta become traceless symmetric products in x-space. These products are 
easily reexpressed in terms of Gegenbauer polynomials. Then the Wick rotated integral 
is solved by introducing spherical coordinates. As a first step one performs all angular 
integrations. Due to the orthogonality of the Gegenbauer polynomials, in this way the 
number of independent summations, stemming from the series expansions, is reduced. 
Finally, the radial integration is performed, and the remaining summations are evaluated 
in order to obtain the final result. 

In particular the GPXT turns out to be very useful in the following situation: 

1. Each angular integration involves maximally two Gegenbauer polynomials. For 
this it might be necessary to expand products of Gegenbauer polynomials with the 
same argument in a Clebsch-Gordan series. 

2. The integrals contain a vertex on which a large number of propagators end, (in our 
case provided by the composite operator) which is chosen as the 'root vertex'. 

3. By partial integrations the momentum dependence in the numerator can be rear- 
ranged such that it is given by a linear combination of scalar products of those 
momenta which run along propagators directly connected to the root vertex. 

4. Only the pole structures of the (logarithmically divergent) integrals are of interest. 

Point [1] is essential for solving the integral. It guarantees that the Racah coefficients 
never appear after having performed the angular integrals. Points[2]to|l]are not essential, 
but significantly simplify the problem. They reduce the number of nested (infinite) sums. 

The first three points are guaranteed if any of the four loops contains at least a 
single propagator which is connected directly to the root vertex and if the non-trivial 
numerators do not generate additional vertices which cannot be reduced by a Clebsch- 
Gordan series expansion. Apart from eventually appearing problems with the non-trivial 
numerators, the requirement holds for all diagrams which become tree level or one loop 
if the root vertex is removed. In our four-loop calculation at critical order, this is 
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guaranteed for all wrapping diagrams and range-five diagrams. At higher loops, it might 
only be spoiled by the increasing complexity of the appearing numerators. Therefore, we 
cannot strictly prove the applicability of GPXT. However, the found numerators even at 
higher loop order suggest that it might hold. In any case, GPXT will in principle allow 
for an analytic determination of the pole parts of whole classes of integrals at arbitrary 
loop order. A practical restriction is only given by the increasing number of domains in 
the radial integrations. 

The fact that we are interested only in the pole structure allows for a dramatical 
simplification. We can neglect the momentum-dependent exponential factor in the x- 
space integrals, i.e. we can set the external momentum of the propagator-type integrals 
to zero. We then have to introduce a cutoff R as upper bound of the integration do- 
main. This regularizes the appearing IR divergences in the integral. Since the poles of 
the logarithmically divergent integrals are independent of the momentum, this proce- 
dure does not alter the poles of the integral from which the sub divergences have been 
subtracted. To ensure this, the IR cutoff has also to be used in the computation of the 
sub divergences. 

GPXT has been proposed in [26] , see also [27, 38] for examples and some important 
relations. We will not review the method here, but only stress an important fact which 
we could not find in the literature. For integrals with more than one derivative on 
a propagator line, special care has to be taken. The corresponding computations are 
briefly summarized in Appendix [Dl We will review the method and the refinements in 
a separate publication [39] . 

4.3 Wrapping dilatation operator 

Having worked out all the integrals presented in Table IC.81 we can now write down 
the wrapping contribution to the four-loop dilatation operator in terms of the chiral 
functions x &s follows 



DJ = -8(2C(3)x(l) - (3C(3) - 5C(5))[x(l, 2) + x(2, 1)] - (1 + 3C(3) - 5C(5))x(l, 3) 
+ (3C(3) - 5C(5)) [x(l, 2, 3) + x(3, 2, 1)] + 3, 2) + ^x{2, 1, 3) 

- ^(1 - C(3))x(2, 4, 1, 3) + - C(3)) [x(l, 4, 3, 2) + x(4, 1, 2, 3)] 




(4.12) 
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Applying it to our operator basis (13.101) we then have 

DJ ^ -8(y + 18C(3) - 30C(5))m . (4.13) 

5 The Konishi anomalous dimension at four loops 

We collect the results from the previous three sections. They can be encoded in a di- 
latation operator at four loops which includes wrapping. It can be obtained by adding 
together the range-four subtracted dilatation operator (13. 9p and the wrapping contribu- 
tion fl4A2l) . We find 

Df^ + DJ = (200 - 16C(3))x(l) - (150 - 24C(3) + 40C(5))[x(l, 2) + x(2, 1)] 
+ (88 + 8e3a + 24C(3) - 40C(5))x(l, 3) 
+ (60 - 24C(3) + 40C(5)) [x(l, 2, 3) + x(3, 2, 1)] 

- (^^ - 2/3 - 4e3, + Ate^b - 2iesc + 4^63^) x(l, 3, 2) ^5^^) 

- (^y - 2/3 - 4e3a - 42636 + 2iesc - 42e3d)x(2, 1, 3) 

+ 4(1 - C(3))x(2, 4, 1, 3) - (10 - 8C(3)) [x(l, 4, 3, 2) + x(4, 1, 2, 3)] 

- (12 + 2/3 + 4e3a)x(2, 1, 3, 2) + (4 - 8C(3))x(4, 1, 3, 2) . 

On the basis (I3.10p it gives 

Df" + DJ ^ (416 - 96C(3) + 240C(5))M . (5.2) 

The non-vanishing eigenvalue is hence found to be 

74 = -2496 + 576C(3) - 1440C(5) . (5.3) 

Restoring the dependence on the coupling constant (II. ip . and including also the contri- 
butions at lower order [28], the final result for the range-four Konishi-descendant up to 
four loops read^ 

7 = 4 + 12A - 48A2 + 336A^ + A^(-2496 + 576C(3) - 1440C(5)) . (5.4) 



^Note that due to a sign change of the x(l) contribution in (I4.12p . the coefficient of the ({3) term 
differs from the one in the first version of [1] . 



Moreover, as remarked in the footnote to (|4.8p . the rational part has been corrected w.r.t to the one 
in the first version of the present paper. It agrees with the one found in [36]. 
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In the literature there also exist several conjectures on the value of this anoma- 
lous dimension [12,13,28] and moreover, as mentioned in the introduction, an explicit 
component field theory computation [29] has been presented following our superspace 
calculation. None of these findings shows an agreement. Therefore it is crucial to have 
an explicit result that can be trusted. 

The aim of this paper has been to present all the steps of our calculation. The super- 
space approach that we have used allowed drastic simplifications: several cancellations 
were automatic and the number of graphs that one really had to compute was reason- 
ably small. With all the ingredients given in this paper an interested reader might want 
to reproduce the whole calculation. 
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A Cancellation of non-maximal range-five supergraphs 



In this appendix we use the result of subsection 12.11 to show that the divergent parts of 
range-five diagrams where one hne interacts with the rest of the graph only through a 
vector line sum up to zero. 

Let us consider a generic three-loop, range-four diagram. We want to extend it to range 
five and four loops by connecting it to a free line from the operator, using a vector 
propagator. There are three possible structures for the part of the three-loop diagram 
which directly interacts with the new line. All the possible diagrams coming from each 
class are shown in figures \A.1\ IA.2I and IA.3I 

• Class A (Figure lAH) : 

As can be seen from the first line of Figure performing part of the D-algebra for 
diagram Sai is enough to obtain the same structure of diagram Sa2, with a different 
sign due to the □ = —p^ cancelling one propagator. Since the two diagrams have 
the same color factor, their divergent parts sum up to zero. 

Now consider diagrams Sas and Sa4- they will produce the same result after D- 
algebra, since in a few steps they can be reduced to the same structure (see lines 2 
and 3 of Figure [Al4|) . However, the two diagrams have opposite color factors and 
so their divergent parts cancel. 
Diagram Sa5 is finite. 

• Class B (Figure 

For diagrams Sbi and Sb2 the discussion is similar to the one for Sas and Sa4- a 
partial Z)-algebra for both diagrams shows that they will produce the same result 
(see lines 1 and 2 of Figure 1X31) . and the opposite color factors make the total 
divergent part vanish. 
Diagram 833 is finite. 

• Class C (Figure lAJj: 

For diagrams Sci and Sc2 we have the same situation as for Sai and Sa2- diagram 
Sci can be reduced to Sc2 through the cancellation of a propagator with a □, which 
gives a minus sign (line 3 of Figure R75|) . This difference in sign is responsible for 
the cancellation of divergent parts, since the two color factors are equal. 
Diagram Scs is finite. 
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(d) Sa4, (e) Sa5 

Figure A.l: Four-loop non- maximal diagrams of class A 




(a) Sbi (b) Sb2 (c) Sb3 

Figure A. 2: Four-loop non-maximal diagrams of class B 
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i?2 






L)2 



(a) 5ci 



(b) Sc2 



(c) 5c3 



Figure A. 3: Four-loop non-maximal diagrams of class C 




Figure A. 4: D-algebra steps showing explicit cancellations for class A 




Figure A. 5: D-algehra, steps showing explicit cancellations for classes B and C 



B Range-five diagrams 



In this appendix we present the results of the computation of range-five diagrams. While 
this calculation would not be strictly required for the subtraction of range-five interac- 
tions from the asymptotic dilatation operator, it has been performed as a check of our 
procedure. Moreover, it allowed us to compute the values of three out of the four e 
coefficients in our scheme. 

The completely chiral four-loop diagrams contributing to the asymptotic dilatation 
operator are shown in Figure IB. 01 They are all range five, except for the last one 
which is range four. Each diagram has a different chiral structure, so the corresponding 
coefficient in the dilatation operator is obtained directly from the value of the diagram. 

The diagrams with structures x(2, 1, 4) (which is the reflection of 3, 4)), 4, 3) 
(reflection of x(l,2,4)) and x(l;4) are shown in flgures IRTI IB. 21 and IB. 31 respectively. 
The results of D-algebra for every diagram in these classes are summarized in tables [KT| 
IB. 21 and IB. 31 cancellations between pairs of supergraphs, as well as diagrams which do 
not produce any divergent part, are explicitly indicated. In all the other cases, a ^ 
indicates that the diagram gives a non- vanishing contribution which must be computed. 
For these relevant diagrams, we also give the relative symmetry factor. In Table IB. 41 all 
these relevant contributions are collected. The integrals Jj are given in Table IB. 51 

Let us deflne C(x(- • •)) the coefficient of structure x(- • •) the dilatation opera- 
tor fl3.6p . According to equation fll.2p . in order to compute C{x{- ■ ■)) from the sum of 
the corresponding diagrams, we have to take the coefficient of the - pole and multiply 
it by -8. We find: 

• C(x(l, 2, 3, 4)) = —10, which is a check 

. C(x(3, 2, 1, 4)) + C(x(l, 4, 3, 2)) = 8/3 -8/3 = 

This must be equal to —2(8 + 263^) from (13.61) . So we obtain e^a = —4 

. C(x(3, 2, 1, 4)) - C(x(l, 4, 3, 2)) = 8/3 - (-8/3) = 16/3 
Since this has to be equal to 4ie3fc, we get e^b = — i4/3 

• C(x(2, 4, 1, 3)) = 2 = 18 + 4e3a which is another check. 

. C(x(2,l,3,2)) = 4-8C(3) 

This must be equal to (—12 — 2/5 — 4e3a), so P = 4^(3) 

• C(x(2, 1, 4)) + C(x(l, 2, 4)) = -4/3 - 20/3 = -8 is a check. 

. C(x(2, 1, 4)) - C(x(l, 2, 4)) = -4/3 - (-20/3) = 16/3 

This allows to compute e^c, since from Sie^t + 4ie3c = 16/3 we obtain e^c = i^/^ 

• C(x(l,4)) = — 4 is a check. 
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Figure B.l: Diagrams with structure x(2, 1,4) 





Fbi * 


1 


Fb4 - 




1 


Fb7 ^ 


-Fb8 






Fb2 * 


1 


Fb5 - 


^ finite 




Fb8 ^ 


—Fbi 






-^53 —Fb6 




Fb6 - 


—Fb3 




Fb9 ^ 


finite 





Table B.l: Results of £)-algebra for diagrams with structure x(2, 1, 4) 




Figure B.2: Diagrams with structure x(l,4,3) 





Fci * 


1 


Fc4 ^ * 


1 


FC7 - 


^ -Fc8 








Fc2 * 


1 


Fc5 — finite 




FC8 - 


. -Fc7 








Fc3 — »■ -Fee 




Fee -Fc3 




FC9 - 


* finite 







Table B.2: Results of D-algebra for diagrams with structure x(l, 4, 3) 




Figure B.3: Diagrams with structure x(l,4) 



Fdi — 


—Fd7 




Fd6 ^ 


—Fd12 




Fdu - 


> finite 




Fd2 ^ 


-Fd3 




Fd7 ^ 


-Fdi 




Fd12 — 


> -Fd6 




Fd3 ^ 


-Fd2 




i^D8 ^ 


finite 




Fd13 — 


> finite 




Fd4 ^ 


* 


1 




— -^DIO 




Fdu - 


> finite 




Fd5 


* 


2 


-fbio — 


—Fd9 











Table B.3: Results of D-algebra for diagrams with structure x(l,4) 



Fai 

Fa2 
Fa3 

Fa4 

Fa5 


^ (<72iV)Vi x(l,2,3,4) 
^ (<72iV)V2 x(3,2,l,4) 
^ (52iV)V3x(l,4,3,2) 
^ (5'iV)V4x(2,4,l,3) 
^ (5'Ar)V5x(2,l,3,2) 


Fbi 

Fb2 

Fbi 


^ (52iV)4(j^ + j2_2J6) x(2,l,4) 
^ -(p2Ar)4j^ ^(2,1,4) 
- \(2.1.4) 


Y^Fb,. ^ -2{g^N)\h xi'^AA) 


Fci 

FC2 
FCA 


^ (52iV)4(Ji + J3-2J7) x(l,2,4) 
^ -(52iV)Vix(l,2,4) 
^ -(52iV)V3x(l,2,4) 


Ei^C**^-2(5'iV)V7x(l,2,4) 


Fda 

Fd5 


^ -2(52iV)4(Ji + J8);^(l,4) 

^ 2(527V)Vix(l,4) 


YFD..^-2{g^N)^JsxiiA) 



Table B.4: Summary of I?-algebra results. 



{Any 



1 1 
+ 
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24£4 24£2 4£. 



J9 = 



(47r)^ 



1 1 



1 



3^3 24^2 3e 



■^3 ^^^^ (4^)8 ( 8£4 + 2£3 Ss^ + Ss) 



J4 



{Any 



5 5 11 



24£4 12£3 24£2 4£ 




J7 = 



(47r)^ 



1 1 
+ 



^^^^^^ (47r)8(l2£2 12s) 




Y- 



(47r)8 V4£2 12£ 



12£4 3£3 12£2 



(47r)8 V 4£y 



Table B.5: Loop integrals for range-five diagrams. The arrows indicate contracted 
spacetime derivatives. 



Four-loop wrapping diagrams with vectors 




Figure C.l: Wrapping diagrams with chiral structure x(l, 2, 3) 





Wbi - 


* 


1 


Wb4 ^ 


* 


1 


Wb7 


finite 






Wb2 - 




1 


Wb5 ^ 


-Wb6 




Wbs 


finite 






Wb3 - 


-> * 


1 


Wb6 ^ 


-Wb5 




Wb9 


finite 





Table C.l: Results of I?-algebra for diagrams with structure x(l, 2, 3) 
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Figure C.2: Wrapping diagrams with chiral structure x(l, 3, 2) 





Wci * 


1 


Wc4 — > finite 






Wc2 * 


2 


Wc5 -Wc3 






Wc3 -Wc5 




Wc6 —>■ finite 





Table C.2: Results of D-algebra for diagrams with structure x(l, 3, 2) 



Dl 




Wi 



D2 




D3 




D4 




D5 




w. 



D6 




Figure C.3: Wrapping diagrams with chiral structure x(2, 1, 3) 



Wdi - 




1 




-> * 


2 


Wd3 - 







Wd4 

Wm 



finite 

-Wn3 
finite 



Table C.3: Results of D-algebra for diagrams with structure x(2, 1, 3) 



El 




Wi 



E2 




E3 




EA 




E5 




Wi 



E6 




w, 



E7 




w, 



E8 




Wi 



E9 




Wi 



ElO 



Weh^ 




We12 = 




w, 



E13 




Wi 



Eli 




Elb 




E1& 




Ell 




w, 



EIH 




E19 




W-, 



E20 




E21 




w, 



E22 




E23 




E24 




W. 



E25 




W, 



E26 




W. 



E27 




Figure C.4: Wrapping diagrams with chiral structure x(2, 1) 



Wei ^ 


-Wew 




We8 




-Wei 




Wew 




-We\2 




We22 


finite 




We2 


* 


1 


We9 




finite 




Wew 




-Weyi 




We23 


finite 




We3 ^ 


-We6 




Wew 




-Wei 




Weu 




-Wew 




We2A 


finite 




We4 ^ 


-We13 




Weu 




* 


1 


Wew 




finite 




We25 


finite 




We5 ^ 


* 


1 


We12 




-Wei^ 




Wew 




finite 




We2Q ^ 


finite 




We6 ^ 


-We3 




We13 




-Wea 




We20 




finite 




We27 


finite 




We7 ^ 


-We8 




Weu 




* 


1 


We21 




finite 











Table C.4: Results of D-algebra for diagrams with structure x(2, 1) 




Figure C.5: Wrapping diagrams with chiral structure x(l,3) 



Wfi ^ 


* 


1 




-Wf% 




Wf2 ^ 


* 


4 


VFf12 ^ 


* 


2 




-Wf5 




W>i3 


-Wfi5 




WFi 


finite 




Wfia 


-Wfi& 




Wf^ 


-Wfz 




Wfi^^ 


-Wf13 




Wf& ^ 


finite 




Wfi% 


-Wfia 




Wf7 


* 


2 


Wfu 


-Wfih 






-Wfu 




Wfis 


-Wfu 






* 


4 


Wf19 


finite 






-Wf21 




Wf2G 


finite 





Wf21 
Wf22 
Wf2Z 
Wf2A 
Wf25 
Wf26 
Wf27 
Wf28 
Wf29 

Wfzo 



-Wfiq 

finite 

finite 

-Wf2.A 

finite 
finite 
finite 
finite 
finite 



Wfzi - 


> finite 




Wf32 - 


* * 


4 


Wf33 - 


* ~Wf25 




Wf3A - 


^ -Wf2A 




Wf35 - 


> finite 




VFf36 - 


* finite 




Wf37 - 


* * 


2 



Table C.5: Results of D-algebra for diagrams with structure x(l, 3) 




Figure C.6: Wrapping diagrams with chiral structure x(l) 




Figure [Q6l Wrapping diagrams with chiral structure x(l) (continued) 



Wgi - 


> -Wg34 




Wg15^ 


finite 




Wg29 - 


> * 


2 


Wg43 - 


> finite 




Wg2 - 


> -Wgio 




Wg16 ^ 


finite 




Wgso - 


> -Wg27 




Wg44 - 


> finite 




Wgz - 


> -Wg5 




Wg17 


finite 




Wg31 - 


> -Wg32 




Wg45 - 


f finite 




Wg^ - 


> finite 




Wgis 


finite 




Wg32 - 


> -TyG31 




Wg46 - 


f finite 




Wg^ - 


. -Wgz 




Wg19 


finite 




Wg33 - 


> finite 




Wg47 - 


> finite 




Wg& - 


> finite 




Wg20 


finite 




Wg34 - 


> -Wgi 




Wg48 - 


> finite 




Wgj - 


^ -Wg25 




Wg21 


finite 




Wg35 - 


> -Wg26 




Wg49 - 


> finite 




Wg% - 


^ -Wg28 




Wg22 


finite 




Wg36 - 


> -Wg38 




Wg50 - 


> finite 




Wg9 - 


> -VFg12 




Wg23 ^ 


finite 




Wg37 - 


> finite 




Wg51 - 


> finite 




Wgio - 


. -Wg2 




Wg24 ^ 


finite 




Wg38 - 


> -Wg36 




Wg52 - 


> finite 




Wgu - 


> * 


2 


Wg25 


-Wgi 




Wg39 - 


> finite 




Wg53 - 


> finite 




WgV2 - 


> -Wg% 




Wg26 ^ 


-Wgss 




Wg40 - 


> finite 




Wg54 - 


> finite 




WGVi - 


> -Wg14 




Wg27 


-VFg30 




Wg41 - 


> finite 










Wgu - 


> -Wg13 




Wg28 ^ 


-Wgs 




Wg42 - 


> finite 











Table C.6: Results of D-algebra for diagrams with structure 



Wax 

Wa2 

Wa^ 




(5^iV)4(74/2)x(2,4,l,3) ^ -2{g^N)^hM 

{g^NYh [x(l, 4, 3, 2) + x(4, 1,2,3)]-. -2{g^N)%M 

{g^N)% x(4, 3, 1, 2) ^ -2{g^N)%M 


E Wa** ^ -2(^2^)4(74 + /2 + h)M 


Wbi 

Wb2 




-ig^N)\h + h + 2/7)[x(l, 2, 3) + x(3, 2, 1)] ^ -2(52Ar)4(73 + + 277)M 
{g^N)%[x{l, 2, 3) + x(3, 2, 1)] ^ 2(32Ar)4/3M 
-(5'A^)'/2[X(1, 2, 3) + x(3, 2, 1)] ^ -2{g^N)%M 

{g^Nfih + h + 2h)[x{l, 2, 3) + x(3, 2, 1)] ^ 2{g^N)\h + h + 2h)M 


E W^B** ^ -2(52iV)4(/7 - /9)[x(l, 2, 3) + x(3, 2, 1)] ^ ~4(52iV)4(lT - /9)M 


Wc2 




-{g^N)\2h + 276)x(l, 3, 2) ^ -2(<72iV)4(272 + 276)M 
2(5'iV)^/2X(l,3,2) ^ 4(52^)472M 


Y.Wc*. ^ -2(52Ar)4/6x(l,3,2) ^ -Aig-'NYhM 


Wdi 

Wd2 




-2((72iV)4(72 + 78)x(2, 1, 3) ^ -4(32iV)4(72 + h)M 
2{g^N)%x{'2, 1,3) ^ Aig-'NfhM 




Weu 
Weu 




-{g^N)\l2 + h + 279) [x(l, 2) + x(2, 1)] ^ 2(^2^)4(72 + 75 + 279)M 
(3^iV)472[x(l,2) +x(2,l)] ^ -2(52iV)472M 

{g^N)\h + h + 2/7) [x(l, 2) + x(2, 1)] - -2(52Ar)4(j3 + 75 + 277)M 
-{g''Nfh[x{l, 2) + x(2, 1)] ^ 2(52iV)473M 


E VF^;** ^ 2(52iV)4(77 - 79)[x(l, 2) + x(2, 1)] ^ -4(^2^)^77 - h)M 


Wfi 

Wf2 

Wfi 
Wf^ 

Wf12 
Wf32 
Wf37 




(<72iV)4(73 + 75 + 477 + 27io)x(l, 3) ^ -2{g^N)\h + h + + 2/io)M 

-2{g''N)\h + h + 2/7)x(l, 3) ^ A{g''Nf{h + h + 2/7)M 

(<72iV)4(74 + 75 + 27n - 47i2)x(l, 3) ^ -2{g''N)\h + h + 2/ii - 47i2)M 

2(52iV)473x(l,3)^-4(^^2^)473M 

-(52iV)473x(l,3)^2(^^2^)473M 

-2(52iV)474x(l,3)^4(52iV)474M 

{g-'Nfhxil,^) ^ -2{g''NfhM 


E Wf» ^ 2{g^N)\lio + 111 - 27i2)x(l,3) ^ -4(32Ar)4(7io + 7ii - 27i2)M 


Wgii 

Wg29 




2{g^N)HM^) ^ 2{g^NYhM 
-2{g^Nfhx{l) - -2(52iV)472M 


EWg** ^ 2{g^N)\h - /2)X(1) ^ 2{g^N)\h - l2)M 



Table C.7: Wrapping contributions at four loops 




Table C.8: Loop integrals for 4-loop wrapping diagrams. The arrows of the same type 
indicate contracted spacetime derivatives 



D Integral calculation via GPXT 

We use the Gegenbauer polynomial x-space technique (GPXT) to compute the four- 
loop propagator-type logarithmic divergent integrals. The integrals are dimensionally 
regularized in D — A — 2e — 2(A-|-1) spacetime dimensions, i.e. the parameter A assumes 
the value 

A = l-£. (D.l) 

We briefly describe our notation. Parentheses around spacetime indices denote a 
traceless symmetric product. For two vectors it is defined as 

xi^x-^ = x^^xl - "^9^^ . (D.2) 

In x-space the integrals can be solved by an expansion of the propagators and traceless 
symmetric products in the numerator in terms of Gegenbauer polynomials as follows 

^ • x,){nr,)i . 

where = xj denote the radial coordinates in Euclidean space, and Xi is the corre- 
sponding unit vector in the direction of Xj. We have also abbreviated 



mm 



= min(rj,rj) , maxj^ = max(ri, r^) . (D.4) 



The Gegenbauer polynomials fulfill 

, r(n + 2A) 

Appearing products of two Gegenbauer polynomials with the same argument have to 
be expanded according to the Clebsch-Gordan series 



i=\m-n\ P_g^ 



2 

i!(i+A)r( '"+^"+' +2A) r( +A)r( '"-^"+' +A)r( '"+^"-' +a) 
r(A)2r(^^^if±^ +A+i)r(i+2Ay r( +i)r( "^-^"+' +i)r(^ri^if^+i) ' 
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where it is important to notice that for positive A the coefficients Dx{m,n,i) never 
become singular. 

For compact notation, we introduce the normahzation factor 



Nx{P,L) 



T{\y 



P-L 



(D.7) 



which depends on the number of propagators P and number of loops L. It arises when 
the momentum integral is transformed via a Fourier transformation into x-space, where 
the integral measure is then rewritten in polar coordinates according to 



1 



d x = -Vtn-ir drdx . 
2 



(D.8) 



Here x is the unit angular vector, and the volume of the D — 1-dimensional unit sphere 
is given by 



27rf 

r(f) 



(D.9) 



D.l Sample calculation of /g with two derivatives 

As an example for the computation of the pole part of integrals with two derivatives we 
show how to obtain /g of Table IC.8[ This integral contains five propagators which do 
not end at the root vertex. 

Shifting the derivatives to the root vertex, the required integral can be rewritten as 




<T>.<J> 




(D.IO) 



We hence have to compute the pole part of the integral 



-A 



2 r(A)'' f A^Xi d^X2 d^Xg d^X4 d^Xs X4 ■ X5 e'^iP-i^3-xi) 



(2% 



X 



4A/ 2 2\1+A(' a2 a2 a2 a2 a2 \\ 
2 l^4^5J l^l2^23^34^45^5lJ 



(D.U) 



After neglecting the exponential function, and introducing an IR cutoff R, the ex- 
pansion of this integral in Gegenbauer polynomials assumes the form 



, oo 

i,j,k,l,r=0 



(D.12) 
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where Rx{i,j, k, I, r) and Axii^j, k, I, r) denote the corresponding contributions from the 
radial and angular integrals. They are given by 



1 1 



Rx{i,j,k,l,r) 



dri dr2 dr^ dr^ dr^ r^r2 r^r^ ^ ^ 
Q (maxi2 niax23 max34 max45 maxsi)'^ 



/ mini2 \ f / min23 \ i / min34 \ 2 / min45 \ 2 / miusi 
Vmaxio/ Vmax9?!/ Vmax-iA/ VmaxAR/ VmaxKi 



V 



.maxi2/ Vmax23/ Vmax34/ Vmax45/ Vmaxsi/ ' 
Ax{iJ,kJ,r)= ^ Dx{lJ,m) dxi . . .dx^C^ixi ■ X2)Cj{x2 ■ xs)C^{xs ■ £4) 



m=\l-l\ 



C^(X4 ■ x^)C^{x^ ■ Xi) 



l + l 

= SijSikSimSir ^ j—T-r^Dx{l,l,i)Ci{l). 

(D.13) 

Evaluating these expressions, substituting them into flD.12p . and expanding the result 
into a power series with negative powers in e, we obtain for the pole part after subtraction 
of the sub divergences 



where the higher negative powers in e indicate that the integral fID.lip contained sub- 
divergences. Substituting the above result into (ID.lOp . we find for the required integral 



(D.15) 



D.2 Sample calculation of In with four derivatives 

As an example for the computation of the pole part of integrals with four derivatives we 
show how to obtain In of Table IC.8I This integral contains six propagators which do 
not end at the root vertex. It is the most complicated one, which we have to compute 
at four loops. 
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Shifting the derivatives to the root vertex, the required integral can be rewritten as 



(D.16) 

We hence have to compute the pole part of the integral 

(D.17) 

where the parameter a has to be introduced to regularize additional divergences which 
arise when we neglect the exponential function and introduce the IR cutoff R. It should 
assume the value a — 0. The parameter a also deforms the (traceless and symmetric) 
products in the numerator as 

3^3 • X(2X2) ■Xe^X2-X3X2-X6- —- , ,. xl X3 ■ Xq . (D.18) 

2[1 — a + A) 

The expansion of the integral is then given by 
A(l + A) 



At{i,j,k,l,r,s) )Rx{i,j,k,l,r,s) , 

(D.19) 



i,j,k,l,r,s=0 

A 



1 - q; + A 
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where the radial and angular integrals are respectively given by 



Rx{i,j,k,l,r,s) = 

Jo 



1 1 



^ dri dr2 dr^ dr4 dr^ drg r^r^r2 ^r^^ ^ rg ^ 
(maxi2 max23 max34 max45 maxsg maxgi)-^ 
mini2 \ 3 / min23 \ f / min34 \ I / min45 \ i / minsg \ 5 / mingi \ f 



maxi2/ \max23/ vmax34/ vmax45/ vmaxsg/ vmaxgi 



Ax{i,j,k,l,r,s) ^5is5mkSmiSmr Y] Y] f. , ,w — rT7^A(l,«,"^)-DA(l,J,m)C4(l) , 

m=|i— 1| m=\j—l\ 

^xihj, k, I, r, s) = 6ij5isSmkSmiSmr . . , .vo/ ^TT^ ^a(1, m)C^{l) . 

^ ' ^ ' m=|i-l| 

(D.20) 

Using the identity 

n i+1 j+1 

I] I] S h{^.3.m) 

'm,i,j=0 m=\i-l\ m=\j-l\ 

n—1 n+1 

= /a(1, 1,0) + ^ 2^ /A(m + a,m + 6,m) + ^ /A(m - l,m - l,m) , 

m=l a,&=±l m=n 

(D.21) 

the integral can be cast into the following form 

I"" = ^^^^7Va(10, 4) (a/a(1, 0) + (A/A(m + 1, m) + A/A(m - 1, m)) 



m=l 



oo 

+ ^{fx{m + l,m- l,m) + fx{m - l,m+ l,m))) , 



m=l 

(D.22) 

where we have abbreviated 

fx{i,j,m) = ——Dx{l,i,m)Dx{l,j,m)C^{l)Rx{i,j,m,m,m,i) , 

+ A)[m + A)^ 

gS.m) = (1 _ ^ ^ ;,)(,\ y^y^^ ^ ^y Dx{l, n^ml)R^i^. ^, m, m, , 
A/A(i, m) = fx{i, i, m) - 5fA(«, m) . 



The explicit forms of the functions which appear in the above expressions read 



^^^(^ (l-a + A)(,n-l + A)4m + A)^ ^"W 

Rx{m — l,m — 1, m, m,m,m — 1) , 

R\{m + l,m + 1, m, m, m, m + 1) , 

A' 

(m-l + A)(m + A)4' 
Rx{m — 1, m + 1, m, m, m, m — 1) 



/A(m - 1, m + 1, m) = — — ^tt/^a(1, m - I, m)Dx{l, m + 1, m)C4(l) 



A^ 

fx{m + 1, m - 1, m) = -—Dx{l, m + 1, m)DA(l, m - 1, m)C^(l) 

(m + 1 + Aj(m + Aj^ 

R\{m + 1, m — 1, m, m, m, m + 1) . 

(D.24) 

It is important to remark that one can set a to its correct value a = in nearly 
all contributions, except in A/a(0, 1), where m = 1. In this case, some integration 
domains of the radial integral R\{0, 0, 1, 1, 1, 0) become divergent as - for a —* 0. This 
divergence is precisely canceled by the prefactor which becomes proportional to a in 
this case. A finite contribution hence remains, and it has to be considered to obtain 
the correct answer for the integral. We will comment further on this in a forthcoming 
pubhcation [39]. 

Expanding the result into a power series with negative powers in e, we obtain for 
the pole part after subtraction of the subdivergences 



Inserting this result into (ID.16P , the result for the required integral is hence found to be 
given by 
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